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1. Introduction 

A many-body system rapidly brought out of some initial state by quenching it to either 
a critical point or into a coexistence region of the phase-diagram where there are at 
least two equivalent equilibrium states undergoes ageing pQ. For ageing systems the 
physical state evolves slowly, non-exponentially and depends on the time since the 
quench was performed and hence time-translation invariance is broken. In addition, 
there holds some kind of dynamical scaling, whether or not the stationary states are 
critical. These aspects of ageing can be conveniently studied through the two-time 
response and correlation functions defined as 

= s- a - 1 h ( t -) ■ h(y) v ^y- XR/z , (i.i) 

h=0 W 

C(t,s) = (0(t,r)0(s,r)) = s- b f c (^) , fc(y) ^ Ac/z , (1.2) 

where the observable 0(t, r) (at time t and location r) is typically taken to be the order- 
parameter (j)(t, r). In this work, we shall also study composite fields such as the energy 
density. We denote by h the field conjugate to O (and when O is the order-parameter 
the conjugate field h is the associated magnetic field). The dynamical scaling forms 
fjl.ipi.2H are expected to hold in the scaling limit where both t,s t m i cro and also 
t s ^> t m j cro , where t m i cro is some microscopic time scale. In writing eqs. (jl.ljll.2j) . 
it is implicitly assumed that the underlying dynamics is such that there is a single 
relevant length-scale L — L(t) ~ t 1//z , where z is the dynamical exponent (the presence 
of another relevant large length scale would break dynamical scaling). Non-equilibrium 
universality classes are distinguished by different values of exponents such as a, b, X c , \r 
(which will depend on the observable O and the field h used and also on whether T < T c 
or T = T c ). For reviews, see [H U 12 U ■ 

In trying to find a systematic approach to determine the scaling functions fn t c it 
has been proposed to generalise the dynamical scaling to a local scale-invariance [61 [7] 
which include the transformation 1 1— > (at + (3)/ ('-ft + 5) in time with a>5 — (3 r y = 1. In a 
field-theoretical setting [H [9] the autoresponse function can be formally rewritten as a 
correlator R(t, s) = (4>(t)cf)(s)) where (j) is the response field associated to 4>. From the 
assumption that both and cf> are so-called quasi-primary scaling operators [TQl Ej (see 
section 3), it follows that 

My) = foy 1+a '- XR/z (y - lr^efc/ - l), (1.3) 

where Q(y) is the Heaviside function which expresses causality, a' is an exponent and 
fo a normalisation constant. A similar explicit, if lengthy, expression can be given for 
the autocorrelation. We refer to [TTI [T2| [13] for recent reviews on the derivation of 
these results and on the numerous examples where these predictions have been tested. 



R(t, 



5(0(t,r)) 



8h(s, r) 
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For our purposes, it is enough to note that most of these tests are done in situations 
where the dynamical exponent z — 2. In particular, almost all existing tests for z 7^ 2 
merely tested the prediction ( 11. 3ft . and this for the order-parameter only, see [HI [12] for 
a detailed discussion. The only exception are a few simple models where z = 4 [TH [15] . 

A fuller picture on the validity of the several technical assumptions which are needed 
for the precise formulation of the theory of local scale-invariance (LSI) can only come 
from more systematic tests of its predictions. To this end, we shall study in this paper the 
ageing behaviour of the spherical model with long-range interactions. It was shown by 
Cannas, Stariolo and Tamarit [33] that for quenches to T < T c , if the exchange couplings 
decay sufficiently slowly with the distance then the dynamical exponent z becomes a 
continuous function of the control parameters of the model and that the scaling forms 
(ll.l|ll.2l) hold for the order-parameter. Here we shall extend these considerations to 
the critical case T = T c and shall further look at the scaling behaviour of composite 
operators (i.e. energy density). Specifically, we shall inquire 

(i) whether dynamical scaling holds, and if so, what are the values of the corresponding 
non-equilibrium exponents ? 

(ii) what is the form of the scaling functions of responses and correlators ? 

(iii) which of the composite operators, if any, transform as quasi-primary fields under 
local scale-invariance ? 

In section 2, we review the exact solution of the kinetic long-range spherical model 
and list our results for the non-equilibrium exponents and the scaling functions for 
the order-parameter and for composite fields. Some of the details are treated in the 
appendix. In section 3, we first show that the presently available formulation [7] of 
local scale-invariance cannot explain our results on the space-time form of the response 
functions when z ^ 2. We then announce some results of a forthcoming paper [17] 
on a general reformulation of local scale-invariance for z / 2 before comparing our 
explicit results with the corresponding predictions of that general theory. In section 4 
we conclude. 

2. Exact solution of the long-range spherical model 

The two-time correlation- and response-functions of the order-parameter in the spherical 
model when quenched either to T = T c or else to T <T C are well-known in the case of 
nearest-neighbour interactions [HI Q3H [201 [ZH [22] . These are also known for the long- 
range model when quenched to T < T c [16J. Here, we shall derive the response and 
correlation functions of the order-parameter and of certain composite operators in the 
long-range mean spherical model quenched to T < T c . 
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2.1. Long-range spherical model 

The long-range spherical model is defined in terms of a real spin variable S(t, x) at time 
t and on the sites cc of a <i-dimensional hypercubic lattice AcZ d , subject to the (mean) 
spherical constraint 



\xeA I 



'.2.1) 



where M is the number of sites of the latticqji The Hamiltonian is given by [2] 

n = ~ J2 J ( x - y^ s * fa - s ^ ' ( 2 - 2 ) 

x,y 

where the sum extends over all pairs (x, y) such that x — y ^ 0. The coupling constant 
J{x) of the model is defined by 

J(x) = (Y!\v\-™) \x\-^\ (2.3) 
VyeA / 

when x ^ and vanishes when x = 0; the summation is over all lattice sites except 
y = 0. The last term in (12.21) . J2 X y J( x ~ v)^xi can & l so absorbed into the Lagrange 
multiplier that imposes the spherical constraint, see below. 

The 'usual' spherical model with short-range interactions is given by J sr {x — y) — 
J X//*(sb) ^y,x+fj,(x), where x + fi(x) runs over all the neighbouring sites of x. When a > 2, 
the relevant large-scale behaviour of the above model, (I2.2p and (I2.3p . is governed by 
this short-range model. Here we shall focus on truly long-range interactions such that 
< a < 2. In this case, the dynamical exponent z = o can be continuously varied 
by tuning this parameter, see [24"t [16] and below. For the equilibrium behaviour of the 
model, consult the classic review by Joyce [24] . 

The dynamics is governed by the Langevin equatio 

5H 



d t S(t,x) 



-i(t)S(t,x)+T}(t,x), (2.4) 

S(t,x) 



where the coupling to the heat bath at temperature T is described by a Gaussian noise 
7] of vanishing average and a variance 

(i](t, x)rj(t' : x')) = 2T5(t - t')6(x - x'). (2.5) 

| For short-ranged interactions, a careful analysis [23] has shown that the long-time behaviour is not 
affected whether (|2.1[) is assumed exactly or on average. 

§ In eq. (|2.4|) . fluctuations in the Lagrange multiplier $(t) are neglected. As pointed out in [22], 
these must be taken into account when treating non-local observables involving spins from the entire 
lattice or if the initial magnetisation is nonzero. Here we are only interested in local quantities and 
use a vanishing initial magnetisation. See |25j for a careful discussion on the applicability of Langevin 
equations in long-ranged systems. 
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The Lagrange multiplier $(t) is fixed by the mean spherical constraint. 

The Langevin equation and the variance of the noise in the Fourier space read 

d t S(t, k) = - (u{k) + S(t, k) + rftt, k), (2.6) 

(2.7) 



(ff(t, k)r}(t', k')) = 2T(27r) d 5(t - t')6(k + k') 



where u(k) = J(0) — J(k). The hatted functions denote the Fourier transform of 
the corresponding functions. In the long- wavelength limit \k\ — > 0, the function 
u(k) — > B\k\ a , where the constant B is given by [16] B = limifei^o (J(0) — J(k))\k\^ cr . 

The solution of the above equation is 



S(t, k) 



-w(k)t 



S(0,k)+ / dre^ T y^T)r](r,k) 



(2.8) 



with the constraint function g(t;T) = exp(2 /^dr^r)). The system is assumed to 
be quenched from far above the critical temperature, hence (S(0,k)) = 0; and the 
spins are assumed to be uncorrelated initially, hence the spherical constraint implies 
(S(0, k)S(0, k)) = (27t) d 5(k + k'). Therefore, the spin-spin correlation function when 
t > s is 



(S(t, k)S(s, k')) = (27r) d 5(k + k')C(t, s; k), 



where 



C(t,s; k) 



,-u(h)(t+s) 



/g(t;T)g(s;T) [ 



1 + 2T [ dTe Mk)T g(r;T) 
Jo 



[2.9) 



(2.10) 



The spherical constraint implies 1 = J A C(t,t;k) and gives g(t;T) as the solution to 
the Volterra integral equation [201 HB] 



(2.11) 



g{t- T) = f(t) + 2T fdrfit - r)g(r; T), 
Jo 

with g(0;T) = 1, and fit) = /(£, 0) is obtained from the function 

f(t;r):= / dk exp (ife • r - 2uo{k)t) , 
where denotes the first Brillouin zone of the lattice A. 



(2.12) 



2.2. Composite operators: Correlations and responses 



We shall now consider not only the spin operator S(t, r) but also some composite 
fields, specifically the spin-squared (spin 2 ) operator and the energy- density operator. 
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We denote the spin and spin 2 operators by 

Ox(t,x) :=S(t,x), (2.13) 
2 (t,x):=S 2 (t,x)-(S 2 (t,x)), (2.14) 

respectively. The energy-density operator is defined as 

O e {t,x):=S{t,x)-(£{t,x)) , 

S{t, x) := J ( x ~ x')S{t, x) (S(t, x') - S[t, x)) . (2.15) 

x' 

These composite operators are defined in such a way that their average value is zero, 
and hence their correlation functions are essentially the connected correlation-functions. 
Also note that since the energy is defined only up to a constant there is no unique 
definition of the energy- density operator. 

The distinction between 2 and O e might be better understood as follows. We 
look into the continuum limit of the energy- density operator, at least for short-range 
model, for we shall later discuss that this operator is not quasi-primary under local 
scale- invariance. In the short-range model, the expression for energy in lattice models 
is usually taken as 

T~t = —J 2J S X S X+ ^ X ), (2.16) 

x,fj.(x) 

where x + n(x) runs over the neighbouring sites of x. In such a case, the energy density 
could be defined as ?(cc) = — JX)u S x S x+tl , which in the continuum limit would reduce 
to ?(#) = — J (2S X + n 2 S x V 2 S x ), where /1 is the lattice constant. But if we had added 
an overall constant E = Af = J2 X then the energy density could be defined as 

e(x) = -jJ2S* (S x+fi - S x ) -> -Jfi 2 S x V 2 S x (1 + OGu)) . (2.17) 

Hence H ST = J2 x e ( x ) = •Jf J ' 2 J2x(^^x) 2 , up to boundary terms. Therefore, for our 
model (12. 2p the two operators C 2 (t,a;) and O e {t, x) must be distinguished. 

The connected two-point correlation functions of the composite operators 

C ab (t, s-x- x') := (O a (t; x)O b (s; x')) (2.18) 

are obtained by making use of Wick's contraction as detailed in the appendix. 
Throughout it is implicitly assumed that t > s unless stated otherwise. As we have 
spatial-translation invariance in our system, we shall find that all two-point quantities 
depend merely on the difference r := x — x' of the spatial coordinates. 

The response functions of the fields {O a (t, x)} to the conjugate fields {h a {t,x)} 



{h}={0} 



(2.19) 
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are obtained by linearly perturbing the Hamiltonion, H — > H — J2 a tx ^a(t, x)O a (t, x), 
as detailed in the appendix. The above defined response function can be interpreted as 
the susceptibility of the expectation value of a field to near-equilibrium fluctuations. 

Finally, we also obtain out-of-equilibrium responses of the fields {O a (t, x)} to local 
temperature fluctuations. This we do by perturbing the noise strength T — > T + 5T(t, x) 
and then evaluating the response functions 



n { p(t,s,;x - x') := . 



5T(s, x') 



(2.20) 

<5T=0 



Let us specify at this point the asymptotic scaling forms that we expect for 
the autocorrelation function C a b(t, s) := C a b(t, s; 0) and the autoresponse functions 
TZ a b{t,s) := TZ a b{t, s; 0) and IZa (t,s) := TZ^\t, s; 0). They are expected to behave 
as 

Cy(t, s) = s~ b »fS(t/s), fS(y) v ~°° y- x ol\ (2.21) 

1Z l3 (t, s) = s-^fHt/s), fi{y) v ~°° y-^' z , (2.22) 

^(t, s) = s-^f^(t/s), fp(y) *K° y-Wf, (2.23) 

in the scaling regime where t, s and t — s are simultaneously large. This also defines the 
nonequilibrium critical exponents a^-, bij, aj , A^, Xq, \r. 

We now write the correlation and response functions of some of the fields {O a (t, x)} 
in terms of the spin-spin correlator C(t, s; r), the constraint function g(t; T) and f(t; r). 
The details of these computations are given in the appendix, while the explicit forms of 
these functions and their asymptotics are spelt out in the next subsection. 

2.2.1. The correlation functions: 

We obtain the following expressions for the non-vanishing correlation functions of the 
composite fields. 

• The spin 2 -spin 2 correlation function is found to be 

C 22 (t, s; r) = (0 2 {t, r)0 2 (s, 0)) = 2 [c(t, s; r)] (2.24) 

For the short-range case, this formula has already been found in [28J. 

• The spin 2 -energy-density correlation functions are 

C 2e (t,s-r) = (O 2 (t,r)O e (s,0)) = ^^-d t [g(t;T)C 22 (t,s-,r)j, (2.25) 

and 

C e2 (t,s;r) =C 2 e(t,s;r). (2.26) 

This is a stronger result than the obvious relation C t2 (t, s; r) = C 2e (s, t\ — r) and 
follows from u>(k) = u(—k). 
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• The energy-density-energy-density correlation function is given by 

Cee(t, s; r) = i^j^d t (g{t; T)C 2e (t, s; r)) . (2.27) 

2.2.2. The response functions: 

For the response functions, we obtain the following expressions. Because of causality, 
in all expressions given below the factor <d(t — s) is implied , where the step function 
Q(t — s) = 1 for t > s, and zero otherwise. 

• Responses to the magnetic field hi(t,x), which are obtained when 7i —>■ 7i — 
J2 tx h{t, x)S(t, x), are given by 



rm - fM f (^ r ) ' (2 - 28) 

TZ 21 {t,s;r) =n tl (t,s;r) = 0. (2.29) 

Responses to the conjugate field h 2 (t,x) of spin 2 operator are obtained when 
7i — > 7i — J2t x ^(t, x)(D 2 (t, x) and are given by 

K 12 (t,s;r) = 0, (2.30) 

K 22 (t, s; r) = 4fc n (t, s; r)C(t, s; r), (2.31) 

n t2 {t, s; r) = -K 22 (t, s; r) d t In / f r) , (2.32) 



The expression for lZ 22 (t,s;r) has already been given in [28] for the short-range 
model. 

Responses to the conjugate field h e (t,x) of energy-density operator are obtained 
when 7i — > TC — J2t x keif, x )@t(t, x) and are given by 

K le (t, s; r) = 0, (2.33) 

n 2e (t, s; r) = 2^rW 9 * T )^ 2 (t, s; r)) , (2.34) 

K*(t, s; r) = ^ff) 9 * [g{t; T)TZ 2e (t, s; r)) . (2.35) 

The spin, the spin 2 and the energy- density responses to temperature fluctuation 
are 

TZ ( p{t,s;r) = 0, (2.36) 

^f ) (t, S ;r) = 2(^ 11 (t, S ;r)) 2 , (2.37) 

K?\t,s-r) = ^-^d t {g{t-T)n 2 {t,s-T)), (2.38) 

respectively. 
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In this section, we first explicitly evaluate in the scaling limit the quantities specified in 
the previous subsection, and then identify the critical exponents and scaling functions. 
The treatment is based on previous results and techniques from [T6], [20] . 

In the late-time limit we can approximate the function u(k) ~ B\k\ a , where 
< cr < 2 [16]. Hence the dynamical exponent in this range of a is given by 



a. 



(2.39) 



Furthermore, the large-time behaviour of f(t) and g(t; T) are as follows. The function 
f(t, x) in this limit becomes 

f(t;x) « B r d/(7 G(\x\t- 1/a ) ; B := f e~ 2m \ (2.40) 

Jk 

Here the scaling function G(|u|t _1 / CT ) for any variable u is defined as 

G(\u\t- l/a ) := B H d/a [ e^" e - 2B l fc l <T ', (2.41) 

Jk 



where J, - ■ ■ — (2tt) d J d d k ■ ■ ■ denotes an integral over ~R d . 
The Laplace transform of f(t) is given by the expression 



f L {p) = -A p- 1+d /° + J2A n (-p) n -\ (2.42) 

71=1 

where the universal constant A = \T(1 — d/o~)\B and the nonuniversal constants 
A n = j Ak (2uj k y n - J k {2B\k\ a )- n , for n = 1,2.... We note that A x = 1/2T C . 

Now the constraint equation (12. lip , upon Laplace transforming, becomes 

^-iMuFY (2 ' 43) 

and is solved in the small-p region using equation (I2.42p . Following a similar analysis 
as done for a = 2 case in [20], we find the large-t limit of the function g(t;T), which is 
given in equations (12.441) . (I2.60p . and (12.761) . This asymptotic constraint function has 
three different forms depending on the quenched temperature and the lattice dimension, 
for a given value of the parameter a. The known case for the short-range model one 
can obtain by taking the limit a — >• 2. The three cases are 

• T < T c : This case was treated in [16] for the spin-spin correlator and the spin 
response. We recover their results and further add other correlation and response 
functions of the composite fields. 

• T = T c , a < d < 2a: To the best of our knowledge the quench to criticality has not 
been treated before. We must further distinguish two critical cases. In the first 
case, d can at most be 4 since a < 2. 
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• T = T c , d > 2a: In this second case of a critical quench, the space dimension d is 
not bounded from above. This case includes the mean-field case. 

We now discuss the large-time behaviour of the correlation and response functions in 
these three cases. 

2.3.1. Case L T <T C 

Since the system exhibits space-translation invariance we take x' = 0. We denote 
y = t/s > 1. The constraint function for T < T c in the large-time limit (16] is 

g(t;T)*B Q (l-pj r d l\ (2.44) 
and hence the spin-spin correlation function for T < T c in the scaling regime reduces to 
C(t, s; k) = (l - B^ 1 s d '° y^° e -m°{t + s)^ (2 45) 

in the Fourier space, or 

C(t, s; r) = C y d ^(y + l)- d/ °G(u), (2.46) 

in the direct space, where C = 2 d ' a (l - T/T c ). Here and below, expressions become 
shorter with the use of the three related scaling variables u, v, and w, where 

u =\r\((t + s)/2)-^ = w(l + s/t)- l '° , 
v = \r\((t- s)/2)- 1/a = w(l-s/ty 1/a , 

w = \r\{t/2)~ 1/a . (2.47) 

The autocorrelation function can now be directly deduced since the scaling function 
G(0) = 1 for r = 0. Hence one reads off, see (12.211) and table [U 

hi=0, = f 1 c \y) = C y d ^(y + ir d ^. (2.48) 

Below we list the remaining expressions in the scaling limit. The autocorrelation 
and autoresponse functions are obtained for the composite operators in a similar way 
as is demonstrated for C(t,s;r) = Cn(t, s;r). The non-equilibrium ageing exponents 
are listed in tables Q] and EJ for future reference. 

We first list the non-vanishing correlation functions. 

• The spin 2 - spin 2 correlator, obtained by substituting equation (12.461) into ( 12.2411 . 
is 

C 22 (t, s; r) = 2Cl y d '\y + \)- 2d ^G\u). (2.49) 
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b 


Ac 


Function 


T <T C 


T = T C 


T <T C 


T = 


T 










a < d <2a 


d>2a 


C n 





d/a-l 


d/2 


3d/2 - a 


d 


C-22 





2d/a - 2 


d 


3d -2a 


2d 


c 2e 


1 


2d/a - 1 


d + a 


3d — a 


2d + a 




2 


2d /a 


d + 2a 


3d 


2d + 2ct 



Table 1. Non-equilibrium exponents b, Ac, as denned in (|2.21[) . for several non- 
equilibrium autocorrelation functions in the long-range spherical model. The exponents 
for the short-range model can be recovered by taking the limit er — > 2. 



Function 


a 




T <T C 


T = T C 


T <T C 


T = T C 


a < d <2a 


d>2a 


K n 


d/a - 1 


d/a-l 


d/2 


3d/2 - a 


d 


n 22 


d/a - 1 


2d/a - 2 


d 


3d -2a 


2d 


n e2 


d/a 


2d /a - 1 


d + a 


3d — a 


2d + a 


n 2e 


d/a 


2d /a - 1 


d + a 


3d — a 


2d + a 


n ee 


d/a + 1 


2d/a 


d + 2a 


3d 


2d + 2a 


n\ 


2d/a- 1 


2d/a - 1 


d 


3d -2a 


2d 


K 


2d/a 


2d /a 


d + a 


3d — a 


2d + a 



Table 2. Nonequilibrium exponents a — a' and as defined in (|2.22p and (|2.23|) . 
for several scaling operators in the long-range spherical model. The exponents for the 
short-range model can be obtained by taking the limit a — > 2. 



• The spin 2 - energy- density correlator, obtained by using equations (12.441 12.49[) in 
(J2J25D, is 

or 2 

C 2e (t, s; r) = —± s- l y dla {y + l)- 1 ""^^ A,G(«), (2-50) 
a 

where, the operator D z is defined as 

D z := zd z + d. (2.51) 

• The energy-density - energy-density correlator, obtained by inserting equations 
(EH I23UD into (E2ZD, is given by 

C ee (t, s;r) = ^ s~ 2 y d/a (y + iy 2 - 2d ^ (D u + d + a) [G{u)D u G{u)\. (2.52) 
a z 

Next we write down the non-vanishing response functions. 

• The spin response function, obtained using equations (12.401 12.441) in (12.281) . is given 
by 

K n (t, s; r) = d s- d ^y d ^(y - iy d/a G(v), (2.53) 
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where C\ = J k exp(— B\k\ a ), and v was denned in eq. (12.471) . 

• The non-vanishing response functions to spin 2 conjugate field, inferred from 
equations (PX3T1 12321) using (EUl 12461 12331) . are given by 

K 22 (t, s; r) = 4C d S - d l°y d l°{y 2 - l)- d '° G(u)G(v), (2.54) 

K £2 (t,s; r) = s-^y^iy 2 - l)- d ^-^-G(u)G(v), (2.55) 

a y-l 

where D v is as given in (I2.5ip and u, v were defined in (12.471) . 

• Responses to the energy-density conjugate field, obtained from equations (I2.34[ 
I2.35P using (I2.44L |2~54|) . are given as follows. 

TZ^t, s; r) = H^L s-^Vtf - 1)^ 
a 

n et (t, s; r) = ^ s-^y d ^(y 2 - 1)^ 

( Dl + aDu 2D U D V D 2 v + aD v \ 

x V¥TTF + + T^IF ) G{u)G{v) - (2 ' 57) 

• The spin 2 and energy-density responses to local temperature fluctuations, obtained 
using equations (12441 12331) in (1237H I23B1) . are 

TZ { p(t, s; r) = 2C\ s- 2d /°y d / a (y - iy 2d ^G 2 (v), (2.58) 

1C 2 

Tlf\t, s; r) = — i s - 1 - 2d ' a y d /' T {y - ly^^G^D^M, (2.59) 
a 

respectively. 
2.3.2. Case Ila: T = T c and a < d < 2a 

For T = T c and a < d < 2a, the constraint function has the form 

g(t;T c ) « (4T 6 M r(-l + rf/a))" 1 r 2+d / CT , (2.60) 
and hence the correlation function in the scaling regime reduces to 

C(t, s; k) = 2T C s y 1 -^ I dz e -B\ k \°(t+s-2s Z ) z -2+d/* ^ ^ 

Jo 

while in direct space is given by 

C(t, s; r) = 2T c Ci y^(y + l)^ £ ^f^^jf^y (2-62) 

n=0 n '\ n ' l°> 

where u is given in (12.471) and the function G n {\\j\t~ 1 / cr ) is defined as 

G n (\t>\r 1/a ) := 4 n r +d/,T J B ( 7 1 f e ik -°e- 2m ° t (B\k\' 7 ) n , (2.63) 

Jk 
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for any variable 0. The spin-response function in this case has the form 

Hn(jt, s; x) = Ci S - d V~ rf/2< % " l)- d/r7 G(v). (2.64) 

To avoid presenting lengthy expressions we write down only the leading behaviour in y 
for the correlators and responses in this case. The spin-spin correlation function in this 
approximation becomes 

C(t, s; r) « 2T c d s^V^^GH, (2.65) 

where T c = T c a/(d — a), and w is as given in ( I2.47p . Setting w and v to zero, we can 
read off the ageing exponents, see tables [1] and [2j 

o u = 6u = --l, A£ = A# = ^ - <r, z = a (2.66) 

The other non-vanishing correlators and responses are given as follows, wherein we 
first list the correlation functions. 

• The spin 2 - spin 2 correlator, obtained from equations ( 12.65L I2T24I) . is given by 

C 22 (t, s; r) « 8T 2 C 2 s 2 - 2d/<T y 2 - 3d/tT G 2 (w). (2.67) 

• For the spin 2 - energy correlator, using (I2.60[ 12.671) in (I2.25p . we obtain 

C 2e (t, s; r) « ^2 s 1 - 2d ^y 1 ' 3d ^G(w)D w G(w). (2.68) 
a 

• Finally the energy- energy correlator, using (I2.60[ 12351) in (12.271) . reads 



4T 2 r 2 

CUt, s; r) « — s-W'y-W'iDv + d + a)[G(w)D w G(w)}. (2.69) 



The non-vanishing response functions are listed below. 

• The responses to the spin 2 conjugate field, obtained using (12.401 [2TM| 12.651) in (\2.31\ 
I2.32p . are given by 

TZ 22 (t, s; r) « 8T c Cf s l - 2d ' a y 2 - 3d ' a G 2 {w), (2.70) 

K e2 (t, s; r) « ^2 8 - 2d ^y 1 - 3d ^G(w)D w G(w). (2.71) 
a 

• The responses to energy-density conjugate field, obtained from ( 12.601 12.701) and 
(I2TMI2T55]) . are 

n 2e (t,s;r) &K e2 (t,s;r), (2.72) 



4T C C 2 



2 



ft ee (M;r) « -^p- S - 1 - 2d/CT 2 /- M/CT ( J D w + c /+a)[G( W )D ?i ,GH].(2.73) 



a 



Lastly, the responses to temperature fluctuations, obtained from (I2.60[ I2.64p and 
(12~371 I2T38D. are 

Tlf \t, s; r) « 2C 2 S - 2d/a y 2 - 3d/(T G 2 (w), (2.74) 
2C 2 

Tlf\t, s-r)^-^- s- 1 - 2d/a y 1 - 3d/a G{w)D w G{w). (2.75) 
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2.3.3. Case lib: T = T c and d > 2a 

For T = T c and d > 2a, the constraint function at large times is 

ff(f;T c )« {4T?A 2 y\ (2.76) 

This is just a constant and does not appear in the correlation and response functions 
to leading order in this large-time limit. In this case, the correlation function in the 
scaling regime reduces to 

c( t , S] k) = ( e -m°its) _ e -B\ kn t +s )} ; (2 77) 

and in the direct space is 

C(t, S ;r) = 2T cCl ~ J^h) , (2.78) 

where G_i is as given in (I2.63[) . 

The spin-response function in this case is given by 

n u (t, s; r) = d s- d ^(y - l)- d ^G{v). (2.79) 

Here again we present only the leading behaviour in y of the correlators and responses. 
The correlation function in this approximation becomes 

C(t, s; r) « 2T C s f(t/2, r) = 2T C C X s^^y^^Giw). (2.80) 

Again we read off the critical exponents after setting v = w = 

a u = &u = - - 1, A^ 1 = A^ 1 = d. (2.81) 
The other non-vanishing correlation functions are given as follows. 

• The spin 2 -spin 2 correlation function, substituting (I2.80P in (I2.24p . is 

C 22 (t, s; r) « 8T 2 C 2 s 2 - 2d/a y- 2d/u G 2 {w). (2.82) 

• The spin 2 -energy correlation function, from (12.761 12.82[ 12.251) . 

C 2e (t, s; r) « s 1 - 2d/ °y- 1 - 2d/,r G(w)D w G(w). (2.83) 

cr 

• The energy-density - energy-density correlation function, from (12.761 [2T831 12.271) . is 



AT 2 C 2 



C et (t,s;r) « S - M / CT 2 /- 2 - M / CT ( J D M , + rf + ( x)[GH^G(«;)].(2.84) 



(7 



The remaining non- vanishing response functions follow. 
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• The responses to spin 2 conjugate field, obtained from (12.401 12.791 12.801) and (12.311 
1232]) . are 

TZ 22 (t, s; r) « 8T c Cf s 1 - 2d/a y- 2d/a G 2 (w) J (2.85) 

TZ e2 (t, s; r) « s- 2d ^y- 1 - 2d ^G(w)D w G(w). (2.86) 

a 

• The responses to energy-density conjugate field, otained from (12.761 [2T85T) and (12.341 
I2.35p . are given by 

n 2£ (t,s;r) &K e2 (t,s;r), (2.87) 

AT C 2 

n ee (t,s;r) « — ^ S - 1 - 2d /^- 2 - M / CT (^+rf+a)[G( W ) J D w G(«;)].(2.88) 

• Finally, the responses to temperature fluctuations, obtained from (12.761 12.791) and 
(12371 |2T38|) . are given as 

TZ 2 (T)(t, s; r) « 2C 2 S - 2d/CT ?T 2d/fJ GV), (2.89) 
of 2 

7^ T) (*, s; r) « — i- s- 1 - 2d/a y- 1 - 2d ' a G{w)D w G{w). (2.90) 

(7 

The exponents of these functions, derived in this section, are collected in tables CD and [2j 
2.3.4- Fluctuation- dissipation ratios 

An important quantity, in particular for the case of critical dynamics, is the fluctuation- 
dissipation ratio of an observable, which is defined as 



X ab (t,s) :=T c n ab (t, S] 0) ^frM;0) ^ 1 (2 Q1) 

and its limit value 

X™ := lim ( lim X ab {t, s)) = lim ( lim X ab (t, s)\ =t/s ) . (2.92) 

For case /, that is for phase-ordering kinetics, it was already known that in the 
quasi-static limit s — > oo but t — s fixed and <C s, the fluctuation-dissipation theorem 
still holds |16j . On the other hand, we obtain in the scaling limit s — > oo and y = t/ s > 1 
fixed that, for all observables considered here 

X n (i, s) = X 22 (t, s) = X 2e (t, s) = X e2 (t, s) = X ee (t, s) = ^^V~ dAr . (2.93) 

CLUq 

For d > o we have therefore in this case that 

Xu = X? 2 = X£ = X™ = X£ = (2.94) 

as expected for a low-temperature phase (recall that for d < a the critical temperature 
is zero [24]). 
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In the case of critical dynamics (case Ha and lib) the limit fluctuation-dissipation 
ratios are universal numbers characterising the critical system [20] . For their calculation, 
we can use directly the scaling limit s — > oo with y = t/s being kept fixed. In case Ha, 
it is convenient to obtain the auto-correlation function C(t,s) by directly integrating 
eq. (12.611) . which leads to 



C(t,s) 



2T c da 
d — a 



s l-d/a yl-d/2a/y 



(2.95) 



Combining this with eq. (12.641) . we get 

X n (t,s)=X n (y) = ~(y + l) 



y-i 



a 



(2.96) 



d - a \2 y + 1, 

Similarly in case lib, using equations (I2.78P and (12.791) . and upon substituting the value 
of GU(0) = aG(0)/(d - a), we find 



X u {t,s)=X n (y) 



y + l 



d/tr s 



(2.97) 



In particular, we see that in the quasi-static limit s — > oo with t — s being kept 
fixed (or alternatively y — > 1), \im y _i Xu(y) — > 1 in both critical cases, such that 
the fluctuation-dissipation theorem holds. Similarly, from the relations (12.24112. 25|2. 271) 
and (12. 31112. 34112. 35p we also have lim^i X 2 2(y) = lim^x X ee (y) = lim. y ^x X 2e (y) = 1. 
On the other hand, and remarkably, the limit fluctuation-dissipation ratio turns out to 
be independent of the choice of the considered observable. We find for y — > oo 



A 22 



VOO 

A 2e 



VOO 

A e2 



XI 



1 — a/d for the case I la 
1/2 for the case lib 



(2.9c 



This reduces to the well-known expressions in the short-range model [20J when z = 
a — > 2. We recall that in [28], a slightly different definition for the energy density was 
used, in which case the value for the corresponding fluctuation-dissipation ratio may be 
different. 



3. Local scale-invariance 



The theory of local scale-invariance (LSI) was developed in a series of papers [HI [TJ [29j 
[HI [30] , using local symmetries to fix the response and correlation functions. For recent 
reviews which focus on different types of applications see [TH [T2J, [13] . For our purposes 
here it is sufficient to just quote a few results. A central concept of LSI are the quasi- 
primary scaling operators [7] , which transform in the simplest possible way under local 
scale-transformations, very much in analogy with the (quasi)primary scaling operators 
of conformal field-theory [lOjijJ A quasi-primary scaling operator <\> is characterised by 

|j Specifically, if X is an infinitesimal generator of a local scale-transformation and <f> a quasi-primary 
scaling operator, 5<j) — —eXifi. Usually, the order-parameter corresponds to a quasi-primary operator, 
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a set of 'quantum numbers' (x,£, pi, /3), where x is the 'scaling dimension' of and p is 
sometimes referred to as the 'mass' of <ft (not to be confused with the lattice constant p 
in section 2.1). 



3.1. Response functions 

For a given dynamical exponent z, LSI yields the following prediction for the response 
function of a quasi-primary operator characterised by the parameters (x,£, p, (5): 

m na i3] 

Irl 



i? L5 '(t, s; r) = ^ 8 p3 R(t, s)J*»® ( 



(t-s) 1 /' 

(,\l+a'-X R /z /, \-l-af 
a) (s ) ' (3 ' 1} 

where the exponents a, a' and Xr are related to the parameters (x,£,/i) via 

1 . . 1 . Ar 2x 2£ 

a + l = - a; + x , a' + l = -(x + 2£ + x + 2£ , — = — + —, (3.2) 

z z z z z 

and the parameters (x, £, /2, /?) characterise the response field 0. The space-time part 
J-^^\p) (where p := |p| and p = r(£ — s) -1 ' 0- ) satisfies the following fractional 
differential equation 

(d p + zpftfT* + W + /i(2 - z)^) F^\p) = 0. (3.3) 

which also illustrates that the 'mass' p may be interpreted as a generalised diffusion 
constant. The fractional derivatives <9" are defined and discussed in [7]. Recall, however, 
that the definition used here is not unique and that different non-equivalent definitions 
for fractional derivatives exist [3U [32]. If z = N + p/q, where N = [z] is the largest 
integer less or equal to z, < p/q < 1 and p and q coprime, the solution of (13. 3p by 
series methods is particularly simple, with the result [1] 



?M (p) = c m <P {m) (p) , with (m) (p) = h n l) p {n ~ 1)z+p /q+m+1 ■ (3.4) 

mg£ n=0 

The constants c m are not determined by LSI and the set £ is 
E _{ -l,0,...,iV-l p^O 

£ -\ 0,...,N-1, p=0 ■ (3 ' 5) 

Finally, the coefficients bn^ read 

b (m) = {-z 2 p) n V{p/q + 1 + m)T(n + z^jp/q + m) + + 2 - z) 
n T((n-l)z + p/q + m + 2)T{z- 1 {p/q + m)+f3 + 2-z) ' { '' 

but if <p is quasi-primary, then neither d t <j) nor d r (j) are. The n-point functions (cf>i . . . (f> n ) of quasi- 
primary operators transform covariantly and hence satisfy linear differential equations X^ (</>i . . . 4>n) = 
0. 
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such that 0( m )(p) has an infinite radius of convergence for z > 1. 

Let us now consider the magnetic response of the order-parameter, TZn, the result 
for which we recall from (12.281) is 

Ku(r;t,s) = (2tt)V^ ( -) F \t/s- l)~ d/ff [ e ife -(*-*r 1/CT & -bv 



R(t,s)J2^p 2n , P = r(t - s)- 1/CT , (3.7) 



n=0 



where the exponent F is given by 

—d/a case (I) 

F = { -2 + d/a case (Ha) . (3.8) 
case (lib) 

Clearly, the space-time part of the LSI-prediction does not agree with this result since 
the exponents of p in eqs. ( 13. 7p and (I3.4p are linearly independent if z is not an integer. 
In eq. ( 13.71) . we have expanded the exponential in order to rewrite this as a series in 
p = \p\. This form of the series is incompatible with the expected form (13.41) for z < 2. 
This disagreement has motivated us to look for a new formulation of LSI, which uses 
a more appropriate form of fractional derivatives V°. This formulation, including the 
exact definition of V" will be described elsewhere in detail [17], here we only mention 
two results we need: 

1. Generalised Bargmann superselection rule: Let a system be given with 
dynamical exponent z ^ §fif> (k £ Nj. Let {&} be a set of quasi-primary scaling 
operators, each characterised by the set (xi,£i, pi, (3i). Then the (2n) -point function 

:= (Mti,r 1 )...<p 2n (t 2n ,r 2n )). (3.9) 

is zero unless the pi form n distinct pairs [pi,p T u\) (i = 1, . . . n), such that 

Pi = -Pr(i)- (3.10) 

2. The decomposition ( 13. lj) of the response function remains valid, but its space-time 
part now satisfies the fractional differential equation, which is quite similar to eq. (13.31) 

(d p + z^pW 2 - z + \f3p + M2 - z)}d p V- z ) F^Xp) = 0. (3.11) 
A solution of equation (13. lip reads [T7J 

r^\p) = f [ e [ <> k |Af exp ( -_l_|fc|* ) (3.12) 
Jk V 2 1 P / 

We see that this prediction of the 'new' formulation of LSI is fully compatible with our 
exact result (13.71) for TZn(t, s; r) if we identify 

p 1 = -p 1 = (z 2 Bi 2 -*)-\ ft=A = 0, g =(2n) d . (3.13) 
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and set for the critical exponents 

ail = a' u = --l, \% = d+^. (3.14) 

This agreement supports the assumption that the fields and are both quasi- 
primary with fi = —p, and (3 = (3. This is further supported by the fact that 
1Zi 2 {t,t';r) = = TZu(t, t'; r), which is predicted by LSI because of the generalised 
Bargmann superselection rule. 

Having verified that the response function for the order-parameter field agrees 
with LSI, and thus having confirmed that is indeed quasi-primary, we now inquire 
whether this holds for composite operators. First, we consider the short-range model 
a > 2. The relevant results can be read from those of section 2 if we let o — > 2. Then 
the response TZn(t, s; r) in eq. (13. 7p simplifies to 

Kn(( , s;r) = (iy 1 (i _ (__L_r!_) , (3,5) 

up to a normalisation constant. Similarly, the temperature response of the spin 2 field, 
from the above expression and eq. (12.371) . becomes 



<W)=*-(|HH~^p(-^), (3.16) 



which is of the form predicted by eq. (13.121) . if we identify 

/i2 = -p-2 = 2/ii, (3 2 = (3 2 = 0, (3.17) 

and 

a 22 = a' 22 = 2a n + 1, A^ 2 = 2A^. (3.18) 

Physically, we can therefore identify temperature changes as the conjugate variable to 
the spin-squared operator, at least for the short-ranged case. On the other hand, the 
spin 2 response 7Z 22 to the perturbation h 2 (t,x) cannot be cast into that form. This 
can easily be seen in equation (12.541) . which has a dependence on t + s, while the LSI- 
predicted form does not contain this dependence. Note that this response function in a 
field-theoretical setting (see for example [HI [9]) corresponds to (0 2 (t, cc)(00)(s, x + r)). 

Our findings suggest that for the short-range model the operator 2 , corresponding 
to spin 2 , is quasi-primary and so is the corresponding response field <fi 2 (obtained by 
locally perturbing the temperature). The parameters of these two fields are related to 
the fields 4> an d <p i n the following way: If <j) has the parameters (x,C,,fi,/3) then the 
parameters of 2 can be obtained from these by multiplying each parameter by the 
factor 2. Similarly the parameters of 2 are related to those of 0. On the other hand, 
we see that the composite operator 00 (defined by a perturbation of the external field 
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h,2(t,x)) is not quasi-primary, and neither is the energy-density operator e(x), even in 
the short-range model (that last finding is not surprising, since we have already seen in 
section 2 that e(x) is related to the gradient of 4>)\% 

We now proceed to the long-range model, where < o < 2. 7^-22 cannot be brought 
into the LSI-predicted form, for the same reason as mentioned above for the short- 

(T) 

range model, namely by comparing the t + s dependence. The response function 7Z 2 
cannot be brought into the LSI-predicted form either, since it contains a product of the 
type T^'{t, s; r) 2 . This again cannot be cast into the general form ( 13.121) . except for 
z = 2. In this exceptional case, the special properties of a Gaussian integral ensure that 
J-^'P^t, s; r) 2 can be rewritten in the form (13.121) upon redefinition of parameters. By 
a similar analysis we find that does not have the LSI-predicted form. We conclude 
that the operator 2 is not quasi-primary under LSI for the long-range model, unlike 
for the short-range case a > 2. 

In a similar way, we also find that the response functions of the operator O e , namely 
lZ t 2 an d 7Z ee , also do not have the form (13.11) and (13.121) . 

Summarising, we have seen that in the long-range model the above composite fields, 
though made of quasi-primary fields, are not quasi-primary. For the time being, the 
order-parameter <fi and the associate response field <ft related to a magnetic perturbation 
remain the only scaling operators with a simple transformation under local scale- 
transformations. This is distinct from the short-range case of z = 2. It remains an 
open question in which sense the transformation of, say, <fi 2 is distinct from the one of 
cf). On the other hand, the generalised Bargman superselection rule (which follows from 
the weaker Galilei-invariance alone) has been confirmed in all cases, by assigning the 
following (relative) 'masses' to the fields 

H = Vo 2 = 2 /i, A*o 6 = 2/x, (3.19) 

and with negative masses to the corresponding response fields. This is natural because 
of the linear structure of the theory. 

3.2. Correlation functions 

In this section, we compare the LSI-prediction for the correlation function of the quasi- 
primary operator (f>(t,x) with our exact result, see (12. 46112. 61112. 77[) . 

The LSI-prediction for the correlation function, for fully disordered initial conditions 
with white noise, is [17]: 

C LSI (t, s; r) = C^(t, s; r) + C^ SI (t, s; r), (3.20) 

In the Landau-Ginzbourg classification of primary scaling operators in the minimal models of 2D 
conformal field-theory (Ising, Potts etc.), one usually has that </> and eventually a finite number of 
normal-ordered powers : 4> : l are primary. 
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with the 'initial' part 

QLSI^j. s - r ^ — c ^ s -b ini t+2f3/z+d/Zy-b init +\ R /z+l3/z fy _ ^b init +d/z-2\ R /z 

x /j fc | 2/3ex p(-i^tl( t + s )) eir ' fc ' ( 3 - 21 ) 
and the 'thermal' part 

(jLSIU^ g . r ^ _ 2T s -hh+2P/z+d/Zy2£/z^y _ ^\2(l+af)-2X R /z-4^/z 

X [ d9 (y - Qy 2 ( a '+ 1 )+ X Rl z + 2 ^./ z +Pl z + d / z [l _ Qy2{a'+l)+\ R /z+2(_/z+f3/z+d/z (3 22) 

X 



Here the function g(w) is not determined by the dynamical symmetries and £ and £ can 
be considered as free parameters. 

In case I, the spin-spin correlation function (12.461) can be rewritten as 

C(t, s; r) = s- F y- F ' 2 [ e ir-* e -BW'(t+.) ) (3.23) 

Jk 

up to a normalisation constant, with a given by (13.81) . In this case (T < T c ), the 
contribution coming from the initial noise is the relevant one [T], and therefore we 
should compare with the spin-spin correlator C^fJ(t, s; r). Indeed we find for the choice 
of parameters as given in (I3.13p . (13.141) and b ini t = that C-£? t (t, s; r) = C(t, s, r), as it 
should be. 

In case Ha and lib, the correlation function, as given in (I2.6ip and (12.771) . can be 
rewritten in direct space as follows, using again (13.81) and up to normalisation constant, 

C th (t,s;r) = 2Tsy- F ' 2 f&6 6 F [ ^1*1^+1-20) e ir-.fc (3 24) 

JO Jk 

For the cases I la (T = T c , a < d < 2a) and lib (T = T c ,d > 2a), in the LSI-prediction 
the term coming from the thermal noise is the relevant one [HE]. If we set g(u) = 1 
and, in addition to the given choice of parameters (13.131) and (13.141) . let 

b th = -- 1, and £ = --zF, f = -zF, (3.25) 
z 4 4 

we find agreement of the LSI-predicted correlation function C t ^ 5/ (t, s; r) = C(t, s; r). 



4. Conclusion 

We have analysed the kinetics of the spherical model with long-range interactions when 
quenched onto or below the critical point T c . For T < T c we have reproduced the 
results of Cannas et al. [16] for the order-parameter and for T = T c we have derived 
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exact results for the response and correlation function of the order parameter. We also 
considered, for T < T c , various composite fields and derived their ageing exponents and 
scaling functions as listed in section 2. We then have carried out a detailed test of local 
scale-invariance using our analytical results. For this purpose, the long-range spherical 
model offers the useful feature that its dynamical exponent z = o depends continuously 
on one of the control parameters. 

We have obtained the following results: 

(i) Dynamical scaling holds for various composite fields for quenches onto or below the 
critical temperature. The non-equilibrium exponents are given in table [1] and [2j 
The scaling functions also have been determined. 

(ii) In the kinetic spherical model with short-ranged interactions (er > 2 and hence 
z = 2), apart from the order-parameter field 0, its square too appears to be a 
quasi-primary scaling operator, as tested through several two-time response and 
correlation functions. 

(iii) In the long-range spherical model, the first tests of the space-time response in a 
system with a tunable dynamical exponent have been performed. This shows that 
the formulation of LSI with z / 2, which we proposed earlier even with the 
recent improvements given in [H], does not describe the exact result for IZn when 
< z < 2, although that formulation did pass previous tests when z = 2 [33] or 

z = a [HUE]. 

(iv) As can be seen from the fractional differential equation satisfied by the space-time 
response function, the precise definition of the fractional derivative used is crucial. 
We shall present elsewhere a systematic construction of new generators of local 
scale-invariance [17] where we shall also show that all previous tests where z = 2 
or z = 4 are passed by the new formulation. Here we have seen that the exact 
results from the long-range spherical model are completely consistent with the new 
formulation of local scale-invariance. 

(v) In contrast to the short-range case where z = 2, the spin-squared field in the long- 
range model is no longer described by a quasi-primary scaling operator. This calls 
for a more systematic analysis, since it indicates that there might be new ways, not 
readily realized in conformal invariance, of non-quasi-primary scaling operators. 

(vi) Both the two-time response and the correlation function of the order-parameter field 
<j) are fully compatible with local scale-invariance in the entire range < z = a < 2. 

While the analytical results presented here certainly provide useful information, the 
eventual confirmation of local scale-invariance might appear fairly natural since the 
underlying Langevin equation is linear. Indeed, for linear Langevin equations there 
is a direct proof of local scale-invariance which uses a decomposition of the Langevin 
equation into a 'deterministic part' for which non-trivial local scale-symmetries can 
be mathematically proven and a 'noise part' [29l [131 Ej- For non-linear Langevin 
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equations the formal proof of non-trivial symmetries of the 'deterministic part' is still 
difficult, although progress has been made [31]. In the absence of exact solutions for 
models described in terms of non-linear Langevin equations numerical tests going beyond 
merely checking the autoresponse function R(t, s; 0) will be required and it will be useful 
to be able to vary the value of the dynamical exponent z. In this context, a natural 
candidate for such studies is the disordered Ising model quenched to T < T c , where it 
is already known that z depends continuously on the disorder and on temperature, see 
|35| [36j [12] and references therein. Furthermore, its Langevin equation is non-linear. 
We hope to be able soon to report tests on the space-time behaviour of response and 
correlators in this model which should provide useful information on whether LSI with 
z 7^ 2 can really be extended beyond the simple solvable systems studied so far. 



Appendix. Correlations and responses 



We present briefly the calculational details that lead to the expressions given in 
section 12.21 In evaluating the expectation value of the composite operators we use 
Wick's contraction, which is applicable in our model if, apart from the noise, the initial 
spin distribution for S(0,x) is also Gaussian [37] . By Wick's contraction and Fourier 
transforming, we get 

C 2e {t, t'; x-x')= ( e Kk+k' Hx - x ') ^ + ^ c ^ t ,. fe)c ,^ t ,, fe / ) 

= ~dt / e^ k+k '>^g(t)C(t, t'; k)C(t, t'; fe'), (Al) 

9(t) A fc ,A fc , 

resulting in equation (12. 25ft . Similarly, 

C ee (t, f ; x - x') = I I e i(*+*'H«--'> (_ Uk + Uw f C (t, t'- k)C(t, t'- k') 

1 ^A fe ,A fe / 

= I e^ h+h '^-^g(t)C(t, t'; k)C(t, t'; k'), (A2) 

gives the equation (I2.27p . 

When H — > 7i — J2 tx hi(t,x)S(t,x), the solution to the corresponding Langevin 
equation is 

3 1 (t,k ] h 1 ) = S(t,k) + -= / dTe^y/tffjhtfak), (A3) 

V9(t) Jo 

where the first term on the right-hand side is the unperturbed solution as given in 
equation (12.81) . Differentiating the above expression with respect to hi(t', x') and Fourier 
transforming back gives TZn(t, t'; x—x'), while both 72-21 (t, t'\ x — x') and TZ e \(t, t'\ x—x') 
vanish since (S(t,k)) = 0. 
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When H — > H — J2t x ^{t, x)S 2 (t, x), the solution to the corresponding Langevin 
equation is 

S 2 (t,k;h 2 ) = S(t,k)+2 [dre^^J^- [ S 2 (r,k';h 2 )h 2 (r,k-k').(A^ 

Jo V 9^ ' ^ A fe' 

Here S 2 {r, k'; h 2 ) on the right-hand side can be replaced by S(r, k'), while evaluating 
the response functions, since the difference is of order 0(/if). Clearly Unit-it'] x — x') 
vanishes since the initial magnetisation is zero. Using the above equation we get, 



'^ h ^m$r)„ = ^^-"'-^(m';^-'^- (A 5) 

and then multiplying it by 2 exp(i(fc + k') ■ x) and integrating over A*, and A fc ' results in 
equation (12.311) . When we multiply equation (IA5I) by 2d; fc exp(i(fe-|-fc / ) -x) and integrate 
over Afc and A fc ', we get 



TZ e2 (t, t'; x-x')=aJ ^rrC(t, t'; x - x') [ u^'^e^*-^, 
V 9(t) JA h 



(A6) 



which is rewritten as in equation (I2.32[) . 

When 7i — > 7i — x h e (t, x)O e (t, x), the solution to the corresponding Langevin 
equation S e (t, k; e) = S(t, k) + 8S e (t, k; e), where 

5S e (t,k;e)= fdre^-'K [ (uj k + w w ) S e (r, k 1 ; e)h e (r, fc-fc').(A7) 

Jo V 9W JA h , 

From the above equation we get 



t=0 



y^ld t (e-^ t - t \-^ +k 'y x ' v W)c(t,t'-,k)) . 



(A8) 



Now multiply this equation by 2exp(z(fc + k') ■ x) and integrate over Afc and A fc / to 
get equation f|2.34j) . If we multiply the equation (1A8[) by a factor {ojk + ujy) then the 
only change in the last term is that d t gets replaced by — d\. Further, multiplying by 
exp(i(k + k') ■ x) and integrating over Afc and A fc / results in equation (j2.35j) . 

When the temperature T — > T + T'(t, x) is shifted then the solution Sx(t, k; T') 
evolves just as given in equation (I2.8p . where the mean is (rf{t, &))t' = 0, but the variance 
becomes 

(rj(t, k)rf(t', k')) T , = (rj(t, k)r}(t', fc')> + 2T' k+y {t)5{t - t'). (A9) 
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This implies (5V(i, k;T')) is independent of T' and 



ST'(t', x') 



8 



(S{t,k)S{t,k')) T , 



T'=Q 




Now multiplying this equation by exp(i(fc + k') ■ x), and then integrating over and 
A fc ' results in equation ( 12. 3 7ft . Similarly, multiplying by exp(i(fc + k ) • x) along with the 
factor {uj k + u k >)/2, and then integrating over A fc and A k > results in equations (12.381) . 
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